Comparison of quantum and classical relaxation in spin dynamics 
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The classical Landau-Lilshitz equation with damping term has been derived from the time evo- 
lution of a quantum mechanical wave function under the assumption of a non-hermitian Hamilton 
operator. Further, the trajectory of a classical spin S has been compared with the expectation value 
of the spin operator S. A good agreement between classical and quantum mechanical trajectories 
can be found for Hamiltonians linear in S respectively S. Quadratic or higher order terms in the 
Hamiltonian result in a disagreement. 
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The Landau-Lifshitz equation [l[ is one of the most of- 
ten used equations in physics. This equation is of impor- 
tance not only in micromagnetism [2] or for the spin dy- 
namics at the atomic level [3| , but also in disciphnes hke 
astronomy Q , biology @ , chemistry , and medicine ■ 
In micromagnetism it describes the motion of a magnetic 
moment in a local magnetic field . The equation of mo- 
tion can be augmented easily by additional interactions 
that can be incorporated into an effective field or, e.g., by 
temperature effects. Moreover, there are many similar, 
equivalent , or alternative approaches namely the Bloch 
equation @], the Ishimori equation 19] and the Landau- 
Lifshitz-Bloch equation fiol] . All these approaches are ca- 
pable to describe magnetization dynamics starting from 
a single atomic spin up to several micrometers. 

Originally, the Landau-Lifshitz equation was intro- 
duced as a pure phenomenological equation yy]. Later 
it has been shown that the precessional term can be de- 
rived by quantum mechanics [ll| , but the damping term 
in the Landau-Lifshitz equation remained phenomenolog- 
ical untill Gilbert proposed to use the Lagrange formal- 
ism with the classical Rayleigh damping instead of the 
original Landau-Lifshitz damping to improve the equa- 
tion, resulting in the Landau-Lifshitz-Gilbert equation 

In this publication I will describe an alternative and 
thereby closing the lack of knowledge: a simple deriva- 
tion of the Landau-Lifshitz equation with damping start- 
ing from the quantum mechanical time evolution will be 
given. Such a derivation provides a deeper understand- 
ing of the underlying mathematics and the connection 
between quantum mechanics and classical physics. 

The derivation starts with the quantum mechanical 
time evolution of the state \4'{t)): 
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Under the assumption of a small time step At we 
can expand the time evolution operator U (t + At, t) ~ 

(i-iHAt/h) +0{At^): 



If we further assume that we have a non-Hermitian 
Hamilton operator: 7^ H we get: 



and the norm: 
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Now we are now looking for a normalized wave function 
and make the ansatz: 
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Then, Eq. ^ can be rewritten 
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Further, with the Taylor expansion: VI — w 1 — |r we 
get: 

m + At))-\m) 
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(8) 

where r is given by Eq. (|4]). In the limit At the dif- 
ferential quotient becomes a differential operator At and 
\ip{t + At)) becomes \ip{t)). Finally, we get the following 
modified time dependent Schrodinger equation: 
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This formula is identical with the equation proposed by 
K. M0lnier et al. [l^ for the calculation of Monte Carlo 
wave functions in quantum optics. 

With H = H - iAf , ?i+ = H + iAf (A e Ej, f hermi- 
tian), and (f ) = {tjj{t)\t\ip{t)) Eq. © becomes: 

i;»A|^(i)) = (H_iA[f-(f)])|^(t)) (10) 

N. Gisin [14] has proposed a similar equation, however, 
with the use of H = H — iAH: 



ih-\m) = iii-m-{m\m) 



(11) 



This special case of Eq. (|T0)) is the quantum mechanical 
counterpart of the Landau-Lifshitz equation as will be 
shown below. 

Eq. ([TT|) can be rewritten as: 



ih 



H-iA 



(12) 



and for the corresponding transposed equation we can 
use the fact that H-'" — H+ = H. Then, the transposed 
commutator is given by 
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(13) 



and therefore the corresponding transposed equation 
d 



(14) 



Now, we are able to write down the corresponding von 
Neumann or quantum Liouville equation [l5[ of the den- 
sity operator p: 



dp 
dt 



dt 



dt 



dt 



p-R 



(15) 



In the Schrodinger picture the time dependence of the 
expectation value (S) is invested in p and in the Heisen- 
berg picture in S, which is in the Schrodinger picture 
time independent: 



ih ^ ^ = ihTi I — S 
dt V dt 



(16) 



and therefore we find under usage of the cyclic change 
under the trace: 



d(S) 
dt 



S;H 



(17) 



Please notice there is still a p included on the right hand 
side of Eq. ([T71) . 

To get the Heisenberg equation we interpret the opera- 
tors in the Heisenberg picture and skip the bra's {ip\ and 



ket's on both sides of the equation ((S) = (V'|S|V'))- 
The expectation values are identical in both pictures and 
therefore we finally get: 
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The problem is, there is still an additional p instead of 
S. For 5 = 1 



2 the density matrix is given by: 
1 



(19) 



The factor ^ is just for the normation because Trp = 1. 
The unity matrix 1 does commutate with H therefore this 
term can be skipped and p is equal to the polarization 
(a) a with the Pauli matrix vector a = (ax, ay,az). Here, 
ari,rj £ {x, y, z} are the Pauli matrices. For general S the 
polarization (a) a has to be replaced by {S)S/S with the 
corresponding spin matrix vector S = (Sx, Sy, Sz) and 

2S 
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Under the assumption of being in a pure state: | (S) | = S, 
and the further assumption that the z-axis is the quan- 
tization axis: (S) — (Sz) — hS we get: 



(21) 



Putting this in Eq. ^TE\\ gives the Heisenberg equation: 

(22) 
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In a previous publication (Tl| I have shown that 

an 
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= S X 
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where the cross product and the gradient directly fol- 
low from the definition of the commutator [l7j and the 
additional term occurs if the Hamilton operator is not 
linear in S„. The double commutator term on the 
right hand side is more complicated. Here, we have 
to know that within the Clifford Algebra S x S = iS 
and (S X S) X H = S X (S X H) holds, which is not 
the case for normal vectors. Here S — (Sx, Sy, Sz) and 
H = (Hx, Hy, Hz) are matrix vectors. Alternatively, we 
can use the following relation of the 5*0(3) Lie algebra: 



X X y = xy = xy - yx [x, y] 



(24) 



where x and y are normal vectors and x and y are 3x3 
skew-symmetric matrices: 



-^3 X2 

xs —xi 

~X2 X\ 



(25) 
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with y and y accordingly. This relation can be proven 
with the aid of the Jacobi identity of the cross product. 
For details and additional depictions see [18]. 

In the limit S* — ^ oo and — > we get the classical 
Landau-Lifshitz equation: 

^ = -— SxHeff + — SX (SxHeff) . (26) 

at /is us 

7 is the gyromagnetic ratio coming from the relation be- 
tween magnetic moment and spin, = ImI comes 
from the normalization, A the damping constant, and 
Hefi — —dH/dS the effective field, with classical Hamil- 
ton function H. 
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FIG. 1: (color online) Magnetization as function of time after 
a gaussian field pulse: comparison of classical (bold lines) and 
quantum mechanical trajectories (thin lines). 
(Dz = 0, ^isB, = 0.1, Tw = 0.02, to = 10, fisB^ = 25.27, 
and A = 0.2) 

To prove the agreement between the TDSE and the 
Landau-Lifshitz Eq. we have performed numerical cal- 
culations. In the following we use a simple single spin 
model. The description is just an example and can be 
extended to systems with N > 1. The corresponding 
Hamilton operator H is given by: 

il^-D,(^S,y -fisBj,^^isB,{t)S^ (27) 

The first term of the Hamiltonian describes a uniaxial 
anisotropy with the 2;-axis as the easy axis. The second 
term represents a static external magnetic field in +z- 
direction. The last term is a time-dependent field pulse 

B,(t) = B-e-K^) (28) 

with gaussian shape to excite the spin. In an experimen- 
tal setup using single atoms such an excitation can be 
realized, e.g., by a current pulse coming from an STM 
(scanning tunneling microscope) tip. In the following we 
investigate the two situations: either (i) Dz = and 
Bz OT vice versa (ii) Dz ^ and Bz = 0. In the 
case (i) all terms of the Hamiltonian are linear in S. In 
the second case (ii) the Hamiltonian contains a quadratic 
term. 
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FIG. 2: (color online) Magnetization as function of time after 
a gaussian field pulse: (a) classical trajectory Srj, (b) quantum 
mechanical expectation values {3^}, r) £ {x,y,z} 
(D, = 0.1, fj,sB, = 0, Tw = 0.02, to = 10, fisB^ = 25.27, 
and A = 0.2) 



In a previous publication we have shown that in 
case (i) under the assumption of a negligible damping 
(A = 0) a good agreement between classical and quan- 
tum spin dynamics can be obtained. Case (ii) shows 
without relaxation a disagreement between classical and 
quantum spin dynamics due to the noncommutativity of 
the quadratic terms. 

Fig. [T] shows the results of the calculation of a single 
spin with relaxation for the case (i) (without anisotropy) . 
Again an excellent agreement between quantum mechan- 
ical expectation values (Sri) rj € {x, y, z} and the classi- 
cal trajectories S*^ of Landau-Lifshitz Eq. is found. In 
case (ii) [Fig. [5] the classical trajectories (a) and quan- 
tum mechanical expectation values (b) disagree. This 
is caused b y th e quadratic Hamilton operator due to the 
anisotropy In this case the commutator [S„; H] leads 
to S X dii/dS plus an additional term of the order of 
h. This additional term vanishes in the classical limit 
S ^ oo and leads to the disagreement between quantum 
mechanical and classical trajectory. In the case of a lin- 
ear Hamiltonian the commutator [S„ ; H] does not lead to 
an additional correction and the both trajectories show a 
perfect agreement. Therefore, we can say that the second 
commutator in the damping term does not produce any 
corrections. 

In the case of the quadratic Haniiltonain with 
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FIG. 3: (color online) Overdamped relaxation of magnetiza- 
tion after a gaussian field pulse. Sa, a & {x, y, z} corresponds 
to classical (bold lines) resp. to quantum mechanical trajec- 
tories (thin lines), (a) Dz — and ^sB^ — 0.1, (b) = 0.1 
and nsBz = 0. {Tw = 0.02, to = 10, fisBS = 25.27, and 
A = 0.2) 



anisotropy we have the additional term 11 1 
iDz/h{S':^, S^^,0) which commutes with S. There- 
fore the deviation between classical and quantum 
trajectory does not come from this term. The question 
is whether this correction comes from the processional 
term — i/?i[S„;H] only or whether the relaxation term 
A/?i[S; [S; H]] also leads to a correction. To answer 
this question and to clarify the effect of the damping 
we compare the trajectories in the overdamped limit 
(A ^ 1). Here we assume that the damping dominates 



the dynamics and skip the processional terms: 
overdamped TDSE is given by: 



d A,^ 



the 



(29) 



and the overdamped Landau-Lifshitz equation by: 
/)S A 

_S X (S X Heff) . 

ot US 



(30) 



Fig. [3] shows the trajectories of the overdamped relax- 
ation process after a field pulse excitation. As expected 
in case (i) = 0, we see a perfect agreement 

between the quantum mechanical and the classical curve. 
In case (ii) 7^ 0, = we find the deviation which 
means that the second commutator also produces a cor- 
rection which modifies the correction which comes from 
the precession term. 

In summary I have shown that it is possible to derive 
the Landau-Lifshitz equation from the quantum mechan- 
ical time evolution of a wave function. This derivation re- 
veals the underlying mathematics and assumptions. Dur- 
ing the derivation we get different presentations in differ- 
ent physical pictures and descriptions. 

In quantum mechanics we can find behavior which can- 
not be described by classical physics like quantum tun- 
neling. In a previous publication |Tl| I have shown that 
quadratic or higher order Hamilton operators do not be- 
have classical, meaning the Ehrenfest theorem does not 
hold in these cases. In this publication this concept has 
been used to proof the damping term. In the case of a 
linear Hamiltonian we see a perfect agreement of classical 
physics and quantum mechanics, but for quadratic and 
higher order Hamiltonians a deviation appears, which 
comes from the damping term. Therefore, the described 
formalism gives us the possibility to compare the classical 
with the quantum spin dynamics. 
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